Branching diffusion in the super-critical regime 
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Abstract 

We investigate the long-time evolution of branching diffusion processes (starting 
with a single particle) in inhomogeneous media. The qualitative behavior of the 
I processes depends on the intensity of the branching. We analyze the super-critical 

Ph ■ case, when the total number of particles growing exponentially with positive prob- 

r~| ■ ability. We study the asymptotics of the number of particles in different regions of 

"j^ i space and describe the growth of the region occupied by the particles. 

2000 Mathematics Subject Classification Numbers: 60J80, 35K10. 

(N : 1 Introduction 
o 

Q ■ The mathematical study of branching processes goes back to the work of Galton and 

Watson [12] who were interested in the probabilities of long-term survival of family names. 
O I Later it was realized that similar mathematical models could be used to describe the 

^ ■ evolution of a variety of biological populations, in genetics [3 El [U [10], and in the study 

of certain chemical and nuclear reactions [131 [11]. The branching processes (in particular, 
branching diffusions) are central in the study of the evolution of various populations such 
as bacteria, cancer cells, carriers of a particular gene, etc., where each member of the 
5^ . population may die or produce offspring independently of the rest. 

In this paper we describe the long-time behavior of the population in different regions 
of space when, in addition to branching, the members of the population move diffusively in 
space and the branching mechanism depends on the location. In particular, we'll consider 
regions Ux(t) centered at a point x{t) = vt which is at a linear (in t) distance from the 
origin at time t. We will be interested in the super-critical case (when the total population 
grows exponentially with positive probability). 

Consider a collection of particles in R'' that move diffusively and independently starting 
with one particle. Besides the diffusive motion, the particles can duplicate with the rate of 
duplication v{x), x G M"^, where x is the position of a given particle and f is a continuous 
non-negative compactly supported function. Both copies start moving independently 
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immediately after the duplication (annihilation of particles and creation of more than 
two particles from one could also be considered, but here we discuss only duplication for 
clarity of exposition). 

Let n^{U) be the number of particles in a domain U C M*^ at time t, assuming that 
at time zero there was a single particle located at x. The large time behavior of n^{U) 
depends crucially on the magnitude of v, that is on whether the operator 

Cu{x) = -Au{x) +v{x)u{x) (1) 

has a positive eigenvalue. This is the operator in the right hand side of the equations 
on the particle density and higher order correlation functions, given below. In the super- 
critical case (i.e., if there exists a positive eigenvalue), if U is fixed, the results of [6] 
on the asymptotics of n^iU) cover, in particular, the case of compactly supported v. 
Namely, n^{U) grows exponentially with a random coefficient in front of the exponent. 
The random coefficients corresponding to different domains differ only by a multiplicative 
constant. (See also |T1], [3]). The distribution of the random coefficient can be described 
in terms of its moments (see [12] )• (For the asymptotic properties of branching random 
walks see also pQ, [2], [IS]-) If has no positive eigenvalues, the total number of 
particles tends to a finite random limit whose distribution can be described in terms of 
its moments. 

In the current paper we consider the super-critical case. We use the spectral techniques 
developed in [5], [12] to get the asymptotic formulas for the density and higher order 
correlation functions of the branching process. These allow us to get the asymptotics of 
fitiUtv), where Utv is a region of fixed size centered at a point whose distance from the 
origin grows linearly with t. This asymptotics is the main result of this paper. 

After recalling the equations on the correlation functions (Section |2]) and the asymp- 
totic behavior of the correlation functions (Section |3]), we show in Section H] that the total 
number of particles, after division by an exponential factor, tends to a random limit in L^. 
In Sections E] and E] we prove a similar result for fixed domains and for domains located 
at a linear (in t) distance from the origin and show that the convergence takes place not 
only in but also almost surely. In Section [7] we show that on the event that the number 
of particles grows exponentially, the region occupied by the particles grows linearly in t. 
In Section IH] we give the distribution of the limiting number of particles in the event that 
the limiting number of particles is finite. 

2 Equations on correlation functions 

Let Bs be a ball of radius 6 in M'^. For t > and x, yi, y2, ... G M'^ with all yi distinct, define 
the particle density pi(t, x, yi) and the higher order correlation functions Pnit, x, yi, ...,yn) 
as the limits of probabilities of finding n distinct particles in Bs{yi),...,Bs{yn), respectively, 
divided by Vor(-B5), under the condition that there is a unique particle at t = located 



2 



at X. We extend p„(t,x, yi, ...^yn) by continuity to allow for yi which are not necessarily 
distinct. For fixed the density satisfies the equation 

dtpi{t,x,yi) = ^Api{t,x,yi) +v{x)pi{t,x,yi), (2) 

pi(0,x,?/i) = 6y^{x). 
Indeed, let s,t > 0. Then we can write 

d f \^ — ^\^ 

+ = (27rs)"2 / e^^pi{t,z,yi)dz + v{x)spi(t,x,yi) + a{s,t,x,yi), (3) 

where the term with the integral on the right hand side is due to the effect of the diffusion 
on the interval [0, s], the second term is due to the probability of branching on [0, s], and 
a is the correction term. The correction term is present since (a) more than one instance 
of branching may occur before time s, and (b) even if a single branching occurs between 
the times and s, then the original particle will be located not at x but at a nearby 
point and the intensity of branching there is slightly different from v{x). It is clear that 
limsio sup ^ y^^d a {s,t,x,y)/s = 0. After subtracting pi(t,x,yi) from both sides of Q, 
dividing by s and taking the limit as s J, 0, we obtain ([2]). 

The equations on p„, n > 1, are somewhat more complicated: 

dtpn{t,x,yi,...,yn) = ^Apn{t,x,yi,...,yn) +v{x) {pn{t,x,yi,...,yn) + Hn{t,x,yi,...,yn)) , 

(4) 

Pn{0,x,yi, ...,yn) = 0. 

Here 

Hn{t,x,yi,...,yn)= ^ piu\it,x,U)pn-\u\{t,x,Y \U), 

where Y = {yi, ?/„), [/ is a proper non-empty subsequence of Y, and \U\ is the number of 
elements in this subsequence. Equation (jl]) is derived similarly to ([2]). The combinatorial 
term Hn appears after taking into account the event that there is a single branching on 
the time interval [0, s], the descendants of the first particle are found at the points in U at 
time s + t, while the descendants of the second particle are found at the points of Y \ U, 
with the summation over all possible choices of U. 

3 Asymptotics of the correlation functions 

First we recall some basic facts about the operator C : L^(]R'^) — > L^(M'^) (see ([T])) and its 
resolvent Rx = {C — A)^^ We will assume that f > is continuous, compactly supported 
and not identically equal to zero. It is well-known that the spectrum of C consists of 
the absolutely continuous part (— oo, 0] and at most a finite number of non-negative 
eigenvalues: 

a(£) = (-CX), 0] U {Xj}, 0<j<N, Xj > 0. 
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We enumerate the eigenvalues in the decreasing order. Thus, if 7^ 0, then Aq = 
maxAj. We assume that there is at least one positive eigenvalue. The resolvent Rx ■ 
L^(M'^) — )■ L^(M'^) is a meromorphic operator valued function on C = C\(— cxd,0]. 

Denote the kernel of Rx by Rx{x,y). If f = (in which case of course there are no 
eigenvalues), the kernel depends on the difference x — y and will intermittently use the 
notations Rx{x, y) and R\{x — y)- The kernel R\{x) can be expressed through the Hankel 
function H^^'': 

Rl{x) = Cdk'^-^{k\x\f-iH^P(iy/2k\x\), k = y/X, ReA; > 0. (5) 

2 

We shall say that / G Cexp(K'^) (or simply Cexp) if / is continuous and 

ll/llaxp(M'^) = sup(|/(x)|el"l') < 00. 

The space of bounded continuous functions on M'^ will be denoted by C(R'^) or simply C. 
The following simple lemma can be found in [5]. 

Lemma 3.1. The operator Rx : Cexp(IR'^) — >■ C{M.'^) is meromorphic in X & C . Its poles 
are of the first order and are located at eigenvalues of the operator C For each e > and 
some A, the operator is uniformly bounded in X E C , |argA| < vr — e, |A| > A. It is of 
order 0(1/|A|) as A— > 00, |argA| < n — e. The eigenvalue Aq of the operator C is simple 
and the corresponding eigenfunction does not change sign. 

From Lemma 13.11 it follows that the residue of Rx at Aq is the integral operator with 
the kernel il){x)il){y), where ip is the positive eigenfunction normalized by the condition 
1 1"*/"! lL2(K<i) = L The function ip decays exponentially at infinity. More precisely, it fol- 
lows from dS]) that if we write x as (6*, in polar coordinates, then there is a positive 
continuous function F such that 

ipi^x) ~ F{6)\x\'2~^ exp(— A/2Ao|a;|) as |a;| — )■ 00. (6) 

For a positive number x, we define the curve T{x) in the complex plane as follows: 

r(x) = {A : |ImA| = V4x(a;-ReA), ReA > 0} |J{A : |ImA| = 2x(l - ReA), ReA < 0}. 

Thus r(x) is a union of a piece of the parabola with the vertex in x that points in the 
direction of the negative real axis and two rays tangent to the parabola at the points 
it intersects the imaginary axis. The choice of the curve is somewhat arbitrary, yet the 
following properties of r(x) will be important: 

First, ReA < x for A G r(x). Second, since the rays form a positive angle with the 
negative real semi-axis, we have |argA| < vr — e{x) for all A G r(a;) for some e{x) > 0. 
Third, since the rays are tangent to the parabola, and the parabola is mapped into the 
line {A : ReA = y/x} by the mapping A — )■ a/A, the image of the curve T{x) under the 
same mapping lies in the half-plane {A : ReA > \/x}. 
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The integration along the vertical lines in the complex plane and along contours T{x), 
below, is performed in the direction of the increasing complex part. 

We'll need estimates on the solutions of the following parabolic equation. Let 

dtp{t, x) = ^Ap{t, x) + v{x)p{t, x), p(0, x) = g{x) G Cexp- (7) 
We'll denote the Laplace transform of a function / by /, 

^ /"OO 

/(A) = / exp(-At)/(t)rft. 
Jo 

Let r be the distance between Aq and the rest of the spectrum of the operator C In the 
arguments that follow we'll use the symbol A to denote constants that may differ from 
line to line. 

Lemma 3.2. For each e G (0,r), the solution of ^ has the form 

p{t,x) = exp{Xot){'ip,g)ip{x) + q{t,x), (8) 

where 

\\qit,-)\\c<AiE)exp{{Xo-s)t)\\g\\c^^^. 
Proof. After the Laplace transform, the equation becomes 

{^A + v)p- Xp = -g. 

Thus, the solution p can be represented as 

p{t,-) = -^ [ e^'RxgdX. (9) 

■^^^ jReA=Ao+l 

The resolvent is meromorphic in the complex plane outside of the interval {—oo, Xq — r], 
with the only (simple) pole at Aq with the principal part of the Laurent expansion being 
the integral operator with the kernel ip{x)ip{y)/{Xo — A). 

By Lemma 13. the norm of Rx does not exceed A/ 1 A | near infinity to the right of 
r(Ao — e). Therefore, the same integral as in (|9]) but along the segment parallel to the real 
axis connecting a point Xo + l + ib with the contour r(Ao — e) tends to zero when 6 — )■ oo. 
Therefore, we can replace the contour of integration in by r(Ao — e). The residue gives 
the main term, while the integral over r(Ao — e) gives the remainder term. □ 

Lemma 3.3. Let K dW^ he a compact set. For each e G (0,r), the function pi{t,x,y) 
satisfies 

pi{t,x,y) = exp{Xot)ip{x)'ilj{y) + qi{t,x,y), 

where 

sup|gi(t,x,y)| < A{e) exp((Ao - e)t - \y\^/2{Xo - e)) (10) 
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for t > 1/2. Moreover, 

qi{t,x,y)ilj{y)dy = (11) 



for each t > 0, x G M.'^, and 

sup / pi{t,x,y)dy < Aexp{Xot). (12) 
Proof. First, let us show that 

(V', Pi{t, y)) = exp{Xot)ilj{y) (13) 
for each t > 0, y E M.'^. Indeed, 

= ^ ((^ + C){exp{-Xos)^),pi)ds = 

(exp(-Aos)^,Pi)|*^o + y ((exp(-Aos)V^),(-— + £)pi)rfs = 

(exp(-Aot)^,Pi(t, ■,?/)) - (V^,pi(0, -,?/)) =exp(-Aot)(V^,pi(t, -,?/)) - V^(l/), 

which proves (IT^ . The relationship (IT^ immediately implies (ITT]) . 

Next, let K' be a compact set that contains supp(f ) U in its interior. In order to 
prove (ITUI) . consider first the case when y E K'. Apply with t replaced by t' = t — 1/2 
and g = pi(l/2, -, y). In order to calculate the main term of the asymptotics, we note 
that lls'llccxp is bounded uniformly in y E K' and 

{ilj,g) = exp(^Ao)V'(2/), 
as follows from ( fT3i) . Therefore, (IHl) implies that 

pi(t, X, y) = exp{Xot)^/j{y)ilj{x) + exp((Ao - e)t)q{t, x, y), 

where \ \qit, ■,y)\\c < ^{K') for all y G K'. We still need to consider the case when y ^ K'. 

Let u{t, X, y) = pi{t, x, y) —po(t, x, y), where po is the fundamental solution of the heat 
equation. Then u satisfies the non-homogeneous version of ([7]) with the right hand side 
/ = —v{x)po(t,x,y) and g = 0. Note that / is a smooth function since y ^ K'. Solving 
this equation for u using the Laplace transform, as in the proof of Lemma 13. 2[ we obtain 



u{t, -, y) = --^ / e^^Rx{-vpo{X, -, y))dX 

-^^^ -/ReA=Ao+l 
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e''Rx{vRU;y))dX (14) 

ReA=Ao+l 



Jr{Xo-e) 



y))d\, 



where the first term on the right hand side is due to the residue at A = Aq. The first term 
can be re-written as 

exp{Xot){tlj,vRl^_^{-,y))ilj{x) = 

exp{Xot){Rl^{v^)){y)ilj{x) = - exp(Aot)V^(i/)^(x). 

The last equahty here follows from the fact that ip is an eigenfunction with eigenvalue Aq, 
that is 

(-A - Ao)^ = -vij. 

In order to estimate the second term on the right hand side of fll4p . we note that from 
it follows that 

\Rl{x,y)\ < A{l)\y/X\^-^x - y\^^'^exp{-V2X\y - x\)\ 
if |A|, |y — x| > /. Thus 

ll^^^(-,2/)lla.,<^(^)ll/|^-^|^l^-^exp(-v/2(Ao-e)|y|) 

for y ^ K', A G r(Ao — due to the fact that Re\/X > VAq — £ for A G r(Ao — and 
\y — x\ > I for X E supp(f), y ^ K'. 

Hence, using the estimate on the norm of Rx : Ccxp — ^ C from Lemma [3. we obtain 

\\R,{vRl{;ymc < A{e)\^\l-I exp{-^2{Xo - e)\y\), A G r(Ao - e). 

Therefore, since ReA < Aq — £ for A G r(Ao — e) and the factor e'^* decays exponentially 
along r(Ao — s), the C-norm of the second term on the right hand side of f|T^ does not 
exceed y4(£:) exp((Ao — e)t — |?/|a/2(Ao — s)). The term po{t,x,y) with x e K, y ^ K', 
t > 1/2, is estimated by the same expression, possibly with a different constant A{e). 
Indeed, if t > 1/2, then 

Po{t,x,y) < Aexp{-\y - x\'^/2t) < Aexp((Ao - e)t - \y ~ x\\/2{Xo - e)) 

since 

\y - x|V2t + (Ao - e)t - \y - xW2{Xo - e) = {\y - x\/V2i - V(Ao - > 0. 

This completes the proof of ( fTOl) . In order to prove ( fT2l) . we again write u{t,x,y) = 
pi{t, X, y) — po(t, X, y). For fixed x, apply the Duhamel formula to the equation 

dtu{t, X, y) = ^Au{t, X, y) + v{y)pi{t, x, y) 
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with the initial data u{0,x, ■) = 0. Now (fT2!) follows since the L^-norm of v{y)pi{t,x,y) 
is bounded by y4exp(Aot) uniformly in x, as follows from (fTOl) . □ 

We'll need additional notations in order to describe the asymptotics of pn with n > 1. 
Let a\{t,y) = ipiri) and al{t,y) = exp{—et — |?/| a/2(Ao — e)). Consider all possible 
sequences a = (ai, cr„) with cXi E {1,2}. By n"(t, yi, we denote the quantity 

U'^{t,yi,...,yn) = sup <H^,2/i)---<"(^,l/n). 

Let Pt : C'cxp — )■ C be the operator that maps the initial function g to the solution 
p(t, ■) of equation ([7]). Let P^g{x) = exp{Xot){ip, g)ip{x) and = Pt — P°. Lemma [3^ 
states that 

IIP/II < A(£)exp((Ao-e)t). 
The particular form of P^P then implies that 

lli^tll < \\Pt\\ + ll^t^ll < ^'exp(Aot). (15) 
For g G C^xp and n > 2, we denote 

In{g) := PnAofi' = / exp{-nXos)Psgds G C. 

Note that ^ ^ 

/ exp{nXQs)Pt-s9ds = exp{nXQt) / exp{—nXos)Psgds 
Jo Jo 

= exp{nXot){In{g) + 0{exp{—{n — l)Xot))) as t — ?■ cxo. (16) 
The functions /i, /2, ••• are defined inductively: fi=ip and 

n-l 



Lemma 3.4. Let K dM!^ he a compact set. For each e G (0,r), the function pn satisfies 
Pnit, X, yi, ?/„) = exp{nXot)fn{x)i){yi) ■ ... ■ ip^yn) + g„(t, x, yi, ...,?/„), (17) 



sup |g„(t,a;, 7/1, ...,?/„) I < y4„(£) exp(nAot)n"(t, yi, ?/„) (18) 



/or t > 1/2. 



8 



Proof. For n = 1, the relation (fT7|) coincides with the statement of Lemma [3.31 Let us 
assume that (ITTI) holds for all natural numbers up to and including n — 1. A generic 
subsequence U CY = {yi, ...,yn) will be written as f/ = (zi, ...,z\u\) and its complement 
as F \ f/ = {zi, ...,Zn-\u\)- By the Duhamel principle applied to the equation for p„, we 
obtain 

p„(t, ^l/l, ...,?/„)= / Pt-s{v p\u\{s,-,Zi,...,Z\u\)pn-\u\{s,-,Zi,...,Zn-\U\))ds 

= / Pt-s{v exp{\U\\os)fiu\i-)i'izi)- ...-^{ziui) 

X exp((n - |[/|)Aos)/„_|;7|(-)V'(^i) ■ ... ■ '4){zn-\u\))ds (19) 
+2/ Pt-s{v ^ exp(|[/|Aos)/|(/|(-)V'(^i) ■ ••• ■^(2:|c/|)gn-|c/|(s,-,^i,---,^n-|c/|))c?s 

+ / Pt-siv q\u\{s,-,Zi,...,z\u\)qn-\u\is,-,Zi,...,Zn-\u\))ds. 

The second and third integrals on the right hand side of (IT^ contribute only to the 
remainder term. Indeed, consider the contribution to the second integral from the term 
with a given U: 

t 

Pt-sivexp{\U\\os)f\u\{-)ip{zi) ■ ... ■ ip{z\u\)qn-\u\{s, ■,zi, ...,Zn-\u\))ds 



< Ail;{z,)...tP{z\u\) / Pt-sivexp{\U\Xos)f\u\{-) 
Jo 

X exp((n - |f/|)Aos)n^-l^l(s, zi, ...,Zn-\u\))ds 

< AiIj{zi) ■ ... ■ ip{z\u\) / exp(Ao(t - s)) exp(nAos)n""l^l(s, zi, ...,Zn-\u\))ds 

Jo 

< Aexp(nAot)^(2;i)...^/'(z|c/|)n^"'^'(t,Zi,...,^„-|c/|) < v4exp(nAot)n^(t, yi, ?/„), 

where the first inequality follows from the inductive assumption and the second one 
from ( fTSl) . The third integral on the right hand side of (fT9ll is estimated similarly. It 
remains to consider the first integral. It is equal to 



..■ip{yn) exp{nXos)Pt^s{v ^ f\u\fn-\u\)ds 
i^iVi) ■ - ■ i'iyn) exp(nAot) (/„(■) + 0(exp(-(n - l)Aot))) 
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where the last equahty follows from (fT6!) . Thus we obtain the main term from the right 
hand side of (fT7|) plus the correction 

ip{yi) ■ ... ■ ip{yn) exp(nAot)0(exp(-(?2 - l)Aot)) 

for which the estimate (fTSi) holds since ip{yi) exp(— Aot) < al(t,yi) due to ([6]). □ 



4 Growth of the total number of particles 

We denote the probability space on which the branching process is defined by (fi, J-", P). 
Let J^t, t > 0, be the filtration generated by the process. We'll write for L^(f2, J-", P). 
Let be the number of particles in M'^ at time t, assuming that at t = there was a 
single particle located at x. 

In this section we prove the basic result on the convergence of Nf /e^°^ in L^. The 
almost sure convergence and the asymptotics of the number of particles in a (possibly 
time-dependent) region of space will be considered in the following sections. 

Theorem 4.1. There is a random variable such that 

a. t^oo, (20) 

where the convergence takes place in LF' . 
Proof. Observe that for < s < t, 

E{N^N^)= / / p2{s,x,yi,z)pi{t- s,z,y2)dzdyidy2+ pi{t,x,y2)dy2. (21) 

jRd jRd jRd jRd 

Indeed, fix yi,y2 G K"^. Then the probability that there is a particle in an infinitesimal 
neighborhood of yi at time s, while a different particle present at time s gives rise to a 
particle in an infinitesimal neighborhood of j/2 at time t is equal to 



P2{s,x,yi,z)pi{t - s,z,y2)dz\ dyidy2. 

/ 

The probability that a particle in an infinitesimal neighborhood of yi at time s gives rise 
to a particle in an infinitesimal neighborhood of y2 at time t is equal to 

pi{s,x,yi)pi{t - s,yi,y2)dyidy2. 

After adding the contributions from the two events and integrating in yi and y2, we 
obtain (EH). 
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Combining ( 12T|1 with ( IT7|) and using that fi = ip, we see that 

E(iVfiVf) = [ [ f [e''"'Mx)ijiy,)ijiz)e'^'^'-Mz)^P{y,) + 

jRd 

6^^°" f2ix)i/j{yi)^p{z)qi{t - s, z, 7/2) + ^2(5, x, yi, z)pi{t - s, z, y2)]dzdyidy2+ 
Pi{t,x,y2)dy2 =: ll/x) + ij^^x) + llt{x) + 



Note that 



since Jj^^ ilj'^{z)dz = 1. Also observe that I]t{x) = since J^^ 4'{z)qi{t — s, z, y2)dz = by 
(fn]). Finally, 

sup \llt{x) + Ilt{x)\ < Ae^o(*+^)-" 
by Lemma [3.41 and ([T^ . Therefore, 

sup |E(iV;iVf ) - e^°(^+*)/2(2;)( / V')'! < Ae^o(*+^)-". 

Thus we have 

\ gAos gAot y g2Aos g2Aot gAo(s+t) 

for X & K. This shows that /e^*^^ is a Cauchy family of random variables as t — )■ cxo, 
and we have convergence in L^. □ 

Remark. It is possible to show (see [I2]) that all the moments of the variables iVf/e^"* 
converge to those of The moments of the limiting distribution are 



ncr=yj^mdyj ux). 

They were shown to determine the distribution of uniquely. 

Remark. In dimensions d > 3 the limiting random variable is equal to zero with 
positive probability. Indeed, since the diffusion is transient, there is a positive probability 
that the original particle wanders off to infinity without branching. Let be the event 
that the number of particles stays bounded (and therefore tends to a finite limit as t — )■ 00) 
and = {^^ > 0} be the event that the number of particles grows exponentially. It is 
possible to show (see [12], for example) that P{B^ U E^) = 1 for each x. 
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5 Growth of the number of particles in a domain 

Let nf{U) denote the number of particles in a domain U C M'^, assuming that at t = 
there was a single particle located at x. Let 



a{U) 



The asymptotics of the number of particles in U is given by the following theorem. 
Theorem 5.1. For each measurable U C ^ we have 

^^^^^a(U)C as t^oo, (22) 

where the convergence takes place in L^. 
Proof. By Theorem 14.11 it is sufficient to prove that 

nf (t/) - a{U)N^ 



gAoi 

Observe that 

TX\2 



— )■ as t — !■ oo. 



E(A^r) = / / P2{t,x,yi,y2)dyidy2+ pi{t, x,yi)dyi, (23) 



E«(f^)^r)= / / P2{t,x,yi,y2)dyidy2+ / pi{t,x,yi)dyi, (24) 



E«(f/)) = / / P2it,x,yi,y2)dyidy2+ / pi{t,x,yi)dyi. (25) 
Ju Ju Ju 

Upon expanding «([/) - a{U)N^y = «(f/))^ - 2a{U)n^{U)N^ + {a{U)N^f, using 
Lemma 13.41 for the asymptotics of pi and p2 and collecting all the lower order terms in 
the remainder R(t,x), we obtain 



E 



h{x){l I i){yi)i){y2)dyidy2 - 2a{U) / / ip{yi)ilj{y2)dyidy2+ 

JuJu JR'iJu 

{a{U)f f [ ij{y,)ij{y2)dy,dy2) + R{t,x) = 

Jw^ jRd- 
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where sup^.^^ x) < Ae This shows that {n^{U) — a{U)N^)/e^°^ is a Cauchy fam- 
ily, thus completing the proof. □ 

Now let us examine the number of particles in the vicinity of a point that is at a linear 
in t distance from the origin. Let b = a/Ao/2 {b will be seen to be the rate of growth 
of the region where the particles can be found with probability that tends to one). Let 
vector with < |f | < 6 and Utv = U + tv the domain obtained from U by 
translation by the vector tv. Let 

e^o* r ih(y)dy 

Note that the asymptotics of g(t) can be obtained from In particular, 

Ait^e(^«'^l''l)* < g{t) < A2t^e(^°-v^l^l)* (26) 

if U is bounded. 

Theorem 5.2. Let U be a bounded domain. For each v eM.'^ such that If I < b, we have 



r as t^ oo, (27) 



9it) 

where the convergence takes place in L^. 
Proof. From (!23|) - (!25|) we obtain 



g{t) e>^otJ ^ ^gi^t) 



Utv 



where 



^(^'^) = (^7a/ / I2{t,x,yi,y2)dyidy2 , ^ . [ [ q2{t,x,yi,y2)dyidy2+ 

Iff If 

q2{t,x,yi,y2)dyidy2 + ^-r- / pi{t,x,yi)dyi- 



Utv 

pi{t,x,yi)dyi + ^— I pi{t,x,yi)dyi). 



9{t)e^^^'Jutv ' e2Aot_. 

The first term on the right hand side of (!28l) is equal to zero as follows from the definition 
of g{t). Applying Lemma 13^ to estimate the integrand in each of the terms in i?(t, x), it is 
easy to see that R{t, x) tends to zero exponentially fast as t — )■ oo. Indeed, let us consider 
the first term. By (1261) and (|T8l) . with an arbitrary e G (0,r), we have the estimate 

q2{t,x,yi,y2)dyidy2\ < 

Utv J Utv 



'9\t) 
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^l-dg2(Ao-v^M)i 

Note that \v\^/2X^ - e - |i;|a/2(Ao - e) < since G (0, a/Ao/2). Therefore, the right 
hand side of (l29l) tends to zero exponentially fast as t — )■ cxd. The other terms in the 
expression for x) can be dealt with in the same fashion. □ 

Remark. With the help of Lemma 13.41 it is possible to show that we have the con- 
vergence of all the moments in f l22|) and f l27j) . 



6 The almost sure convergence 



Theorem 6.1. The convergence in [2D^) takes place almost surely. If U is a domain with 
a smooth boundary, then the convergence in ( fl^) and (21) takes place almost surely. 



Proof. We will only prove the almost sure convergence in fl271) since the other statements 
can be proved similarly. Fix an arbitrary 5 > and K > Q. By the Borel-Cantelli lemma, 
it is sufficient to demonstrate that there is an increasing sequence t„ — )■ oo such that 

f;p( sup i^^i^-rixJ, r<i^)<oo. 

„=i te[t„,i„+i] 9V') 

From the proof of Theorem 15.21 it follows that n^{Utv)/g{t) converges to in expo- 
nentially fast. Let 7 > be such that 

E(^:^i^-n'<^e-^* (30) 
9{t) 

for some constant A. We take t„ = 21nn/7, n > 1. By the Chebyshev inequality, 

^' git^) ^ ' 5^ - 52 52^2' 

and therefore 
It remains to show that 



> P sup 77^ > ^ <is: <oo. 
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which is equivalent to the following two inequalities holding at the same time 

> P sup — — " > ^ <K)<oo, 32 

> P sup 77^+ s > ^ <K)<oo. 33 

We will only prove ( 132|) since ( l33l) can be proved similarly. Let us show that the term 
''^tni^tnv)/ g(tn) in (132!) can be replaced by a more convenient expression. For r > 0, let 
= {x G M'^ : dist(x, f/) < r} be the r-neighborhood of U. Let (/^(t) = g{U^,t,v). 
Since ?7 is a smooth domain, from the definition of g it follows that r can be chosen to 
be sufficiently small so that 

| «''"')-«M |<V^) (34) 

for all sufficiently large n. Moreover, since r > and t„+i — 1„ — )■ oo, we have 

U f^*'' ^ (35) 

tG[tn,tn + l] 



for all sufficiently large n. As in (IHT]) . we have 

ft/'" 



EP(l^^^^-ri>^)<oo. (36) 



n=l 



Now, 



yp( sup |!!l4:i-^_!!^^|>-, e^<ir)< 

a(t) q(U ' 5' ' 



n=l 
oo 



git) 


gitn) 




^liUt^v) 


g'itn) 


gitn) 


1 ^ ' 





The first series on the right hand side is finite by (1311) and (1361) . The second series is 
estimated from above by 

Hn: sup iK+'-)f^^-i^^>h^tn\"^^-a>h■ 

■ ' ^ git) 5 ^ g'{tr^ ^ 



n=l 
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The first term of this expression is finite as follows from (IMIl . while the second one is 
finite by (!36|) . Therefore, (132|) will follow if we demonstrate that 



> P( sup j—^ ^ > -, ^ <K)<oo. (37) 



Similarly to ( 13T]) . we have 



n=l 



for each constant c > 0. Combining this with f l26|) and fl35|) . we see that fl37|) will follow 
if we show that for each a, i? > we have 



^P( sup - nliUlJ > n^ Nl < i?e^"*") < oo. (38) 

n=l *S[tn,in+l] 

Roughly speaking, we need to show that the number of particles that visit a smaller 
region U^^^ over a short interval of time [tn,^n+i] can't significantly exceed the number of 
particles that are in the larger region U^^^ at the time Let us fix n and study the n-th 
term in the series (138|) . After conditioning on J-^^, the questions becomes the following: 
Suppose we have m < Re^°^" particles at time zero located at xi,...,Xm- Let yi,...,yM 
be their descendants at time t = tn+i — tn- Each of these has a starting point in the 
set {xi, ...,Xm}- We are interested in the probability that at least n" descendants of the 
original particles were at a distance r/2 from their respective starting points prior to the 
time t = tn+i - tn- 

The expected number of descendants of a single particle that cover distance r/2 (from 
the initial position of the particle) in time t decays faster than e~^/^ as t — )■ for some 
/3 > 0. Therefore, by the Chebyshev inequality, the n-th term in ( 138|) is estimated from 
above by Re^otnQ-P/itn+i-tn)j^-a^ which decays exponentially in n, as follows from the def- 
inition of tn- Therefore the series (138|1 converges, which completes the proof. □ 



7 Limiting shape of the region occupied by particles 

Let B{r) denote the ball of radius r centered at the origin. Recall that b = 

Theorem 7.1. For each 6 > 0, there exists a random variable T = T{6) (T < oo almost 
surely) with the following properties: 

(a) There are no particles outside B{{b + 6)t) for t >T. 

(b) On the event > the union of the unit neighborhoods of the particles cover 
B{{b - 5)t) for all t > T. 
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Sketch of the proof. We'll only verify the statement for a sequence of times tn = c\nn for 
a certain c. Namely, we'll show that there is a random variable and a constant c > 
such that: 

(a') There are no particles outside of B{{b + 5)t„) for n > N. 

{b') On the event > the union of the unit neighborhoods of the particles cover 
B{{b - 5)t„) for n>N. 

The transition from the sequence of times to the continuous time can be then accom- 
plished similarly to the way it was done in the proof of Theorem 16.11 

By the Chebyshev inequality, 



F{nl{B{{b + 6)tn)) > 1) < Enl{B{{b + S)t^)) = / p,{tn,x,y)dy < 

JB({b+5)t„) 



Ai / {exp(Xotn-\y\V2^) + exp{{Xo - e)tn - \y\^/2jxi^^^))dy < 

JB({h+5)tn) 

A2 (exp(t„(Ao - v^(& + ^))) + exp(t„(Ao - e - [b + 5)^2{\^ - e)))) , 



where the second inequality is due to Lemma 13731 and (j6]). By choosing a sufficiently small 
e > 0, we can make the right hand side of the last formula smaller than A2 exp(— t„5-\/Ao). 
Therefore, 

00 00 00 

n<SB{{b + 5)tn)) > 1) < 5^ e-*"^^ < 5^ < 00 

n=l n=l n=l 

if we choose t„ = cilnra with ci > 2/{5\/\q). By the Borel-Cantelli lemma, there is a 
random variable A^^i such that (a') holds (with A^^i instead of N). 

In order to establish (6'), note that for each n, the ball B{{b — S)tn) can be covered 
by the balls B^, i?^*-"^ of radius 1/2 centered at t„f ^, t„f™*^"^ in such a way that the 
centers of the balls are inside B{{b — 6)tn) and ■m{n) = 0{t'^). Let = g{B{l/2),tn,v''). 
As in (|30|) . there is 7 > such that 

Tn 

By the Chebyshev inequality, for each A; = 1, ...,m{n) and each > 0, 
P«(5^) = 0, C>K)<F{nUB':i)<Kg'j2, C > K) < 



7^2 
Therefore, 

00 00 
^P«(S^) = for some fc, > K) <Y,MK)e-'^'<- 

n=l n=l 

The series on the right hand side of this inequality converges if we choose t„ = C2 In tt, 
with C2 > 2/7. Since K > Q was arbitrary, by the Borel-Cantelli lemma, (6') holds (with 
N2 instead of N). It remains to take c = max(ci, C2) and then N = max(A^i, N2). □ 
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8 Limiting distribution in the case of finitely many 
particles 

In this section we make a couple of remarks concerning the distribution of the total num- 
ber of particles on the event = 0. 

Remark. Let = lim(_>oo-^f- This random variable is finite almost surely on the 
event = (this can be proved similarly to the way it was done for the corresponding 
statement in the critical case in [12]). In other words, with probability one, the total 
number of particles either grows exponentially or tends to a finite limit. The latter event 
has nonzero probability if and only if o? > 3. 

Remark. Let d > 3 and define M"(x) = P(iV^ = n), n > 1. The quantities M"(x) 
satisfy a recursive system of partial differential equations. Namely, 

^AM\x) = v{x), (39) 

with the condition at infinity 

lim M^{x) = 1. 

\x\—^oo 

For 72 > 2, we have 

^ n—l 

-AM"(x) = v{x) M\x)M''-\x), (40) 

k=l 

lim M'"(x) = 0. 

\x\—^oo 

Equations fl39|) and fj40l) can be obtained by considering the behavior of the initial particle 
on the time interval [0, 6] such that 6^0, with the left hand side accounting for the 
diffusive motion and the right hand side for the branching. 
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